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ABSTRACT 

We present the results of an observational study of the efficiency of deep mixing in globular cluster 
red giants as a function of stellar metallicity. We determine [C /Fe] abundances based on low-resolution 
spectra taken with the Kast spectrograph on the 3ni Shane telescope at Lick Observatory. Spectra 
centered on the 4300 A CH absorption band were taken for 42 bright red giants in 11 Galactic globular 
clusters ranging in metallicity from M92 ([Fe/II]= -2.29) to NGC 6712 ([Fe/H]= -1.01). Carbon 
abundances were derived by comparing values of the CH bandstrength index S2{CH) measured from 
the data with values measured from a large grid of SSG synthetic spectra. Present-day abundances are 
combined with theoretical calculations of the time since the onset of mixing, which is also a function 
of stellar metallicity, to calculate the carbon depletion rate across our metallicity range. We find that 
the carbon depletion rate is twice as high at a metallicity of [Fe/II]= —2.3 than at [Fe/II]= —1.3, 
which is a result qualitatively predicted by some theoretical explanations of the deep mixing process. 
Subject headings: Globular clusters: individual (NGC 4147, NGC 5727, M3, NGC 5904, M5, NGC 
6205, M13, NGC 6254, MID, NGC 6341, M92, NGC 6535, NGC 6712, NGC 6779, 
M56, NGC 7078, M15, NGC 7089, M2) - Stars: abundances - Stars: evolution 
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1. INTRODUCTION 

The observation that carbon abundance in globular 
cluster red giants declines continuously as the stars 
evolve has inspired a great deal of observational and 
theoretical study (e.g., Suntzeff 1981, Carbon et al. 
1982, Trefzger et al. 1983, Suntzeff & Smith 1991, 
Weiss & Charbonnel 2004, Denissenkov & Tout 2000, 
Smith & Briley 2006, and similar work). Canonically, 
abundances should be static on the red giant branch 
after the first dredge-up because of the broad radia- 
tive zone between the hydrogen-burning shell and the 
surface. Progressive carbon depletion with rising lu- 
minosity on the giant branch is commonly interpreted 
as a sign of a non-convective "deep mixing" process 
that mixes carbon-depleted material from the hydrogen- 
burning shell, where the CN(0) cycle is acting, to the 
surface (e.g., Sweigart & Mengel 1979, Charbonnel 1995, 
Charbonnel et al. 1998, Bellman et al. 2001, Denis- 
senkov & VandenBerg 2003, and others). This same de- 
pletion of surface carbon abundance is also observed in 
red giants in the halo field (e.g., Gratton et al. 2000), 
and is observed to occur at the same rate in the halo 
field as in globular clusters with halo-like metallicities 
(e.g.. Smith & MarteU 2003). 

The process of deep mixing, inferred from observations 
of low [C/Fe], loge(Li), and ^^C/^^C, is only observed to 
occur in stars brighter than th e red giant branch (RGB ) 
luminosity function "bump" (jCharbonnel et all Il998l ) . 
There are indications (e.g., Langer et al. 1986, Bellman 
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et al. 2001) that there may be carbon depletion in stars 
fainter than the RGB bump in the metal-poor globu- 
lar cluster M92, though there is not an obvious physical 
explanation for such a phenomenon. During the first 
dredge-up, in the subgiant phase, the base of the convec- 
tive envelope drops inward to smaller radius as the stellar 
core contracts, mixing the partially-processed material of 
the stellar interior with the unprocessed material at the 
surface. As hydrogen shell burning progresses, low on 
the giant branch, the temperature gradient in the star 
steepens and the base of the convective envelope begins 
to move outward, leaving behind a sharp jump in mean 
molecular weight (the "/i-barrier" ) at the point of its fur- 
thest inward progress (|Ibenlll965[ ) . As the convective en- 
velope retreats outward, this steep /i gradient finds itself 
within a radiative region between the hydrogen-burning 
shell and the base of the convective zone, where it can po- 
tentially hinder mass motions within the radiative zone. 

The red giant branch bump is an evolutionary stut- 
ter that occurs when the hydrogen-burning shell, which 
is advancing outward in mass, encounters the /i bar- 
rier. The sudden influx of hydrogen-rich material to the 
hydrogen-burning shell causes the star to become briefly 
bluer and fainter before it re-equilibrates and continues 
to evolve along the red giant branch (e.g., Iben 1968, 
Cassisi et al. 2002). In a collection of stars with equal 
age and composition, this evolutionary loop will result 
in an unexpectedly large number of stars at a particu- 
lar magnitude, and a bump in the differential luminos- 
ity function. At a fixed mass, the base of the convec- 
tive envelope sinks lower in higher-metallicity stars dur- 
ing the first dredge-up, meaning that the RGB bump 
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occurs at a fainter luminosity on the RGB in high- 
metallicity globular clusters than in low-metallicity clus- 
ters (e.g., Zoccali 1999). However, because evolutionary 
timescales shorten as stellar mass rises, there is a max- 
imum mass of ~ 2Mq for stars to experience this evo- 
lutionary loop: above that mass, the hydrogen-burning 
shell does not move outward far enough to cros s the iJ. 
barrier in the short time the star is on the RGB (iGilrovl 
[1989). The fact that deep mixing does not begin un- 
til after the hydrogen-bur ning shell crosses t he fi bar- 
rier is interpreted by, e.g., ICharbonnell (|1994| ) to mean 
that the gradient of mean molecular weight is the domi- 
nant factor in permitting or prohib i ting d eep mixing. In- 
deed, iDenissenkov fc VandenBer^ (|2003D point out that 
V/x is "the only physical quantity that changes signif- 
icantly while approaching the hydrogen-bu r ning shell" 
in post-bump red giants. IChaname et al.1 (|2005f ) pro- 
vide a somewhat different perspective: in their maximal- 
mixing models, the fj, gradient inhibits mixing on the 
upper giant branch, but rotational mixing processes are 
not strong enough on the lower giant branch to cause 
observable changes in surface abundances, regardless of 
whether there is a steep gradient present. 

The underlying physical reason for deep mixing is not 
clear , though r otation has commonly been implicated 
since ISweigart & : Mongol (1979) proposed meridional cir- 
culation as an explanation for CNO anomalies in red 
giants. Recent theoretical studies tend to focus on spe- 
cific parametrizations and representation s of the process; 
for example, Rayleigh- Taylor instability (Egglcton et al. 
[2008), diffusion (Dcnisscnkov fc VandcnBcrg 2003) , 
and thermoha l ine m ixing (jCharbonnel fc ZahnI l2007f ) . 
iPalacios et al.l (|2006l ) demonstrated that meridional cir- 
culation, differential rotation, and shear turbulence do 
not create enough mixing to account for the observed 
variations in surface abundances, implying that addi- 
tional hydrodynamical processes must be acting. The 
large stu dy of surface a bund ances in field giants pub- 
lished bv lGratton et al.l (|2000( ) is a key to differentiating 
between models of deep mixing, since it demonstrates 
clearly the progressive depletion of carbon on the giant 
branch, as well as the sharp drop in ^^C/^'^C and log e(Li) 
that happens at the RGB bump. 

The fundamental result from iGratton et al.l (|2000l ) 
that all current deep mixing models must reproduce is 
that deep mixing is universa l among post-bump red gi- 
ants. ICharbonnel et aP (|1998 ) consider mixing in terms 
of the "critical /x gradient," the largest gradient in mean 
molecular weight that still permits deep mixing. In an 
observational study of seven mildly metal-poor red gi- 
ants in the region of the RGB bump, they find that the 
c ritical /x gradie nt is independe nt of c omposition or mass. 

IDenissenkov fc VandcnBer^ (|2003f ) use the formalism 
of diffusion, with mixing depth and a diffusion con- 
stant as the important parameters, to model deep mix- 
ing. They find that the mixing depth does not depend 
strongly on metallicity, which implies that all red giants 
with a mass less than ~ 2M0 will experience deep mix- 
ing. Their figures also show that the evolution of surface 
abundances of carbon and nitrogen are not particularly 
affected by metallicity, thou gh reductions in l oR ^( Li) and 
12(-i^13q g^j.g more sensitive. Eg gleton et al.l ([2008) show 
that the reaction ^He('^He, 2p)'*He causes a ^-inversion in 
the outer edge of the hydrogen-burning shell, and claim 



that the resulting Rayleigh- Taylor instability is i mpor - 
tant in driving deep mixing. ICharbonnel fc ZahnI ([2003) 
argue that the more complex process of thermohaline 
con vection will a ct in that /x-inversion region. They use 
the lUlrichI (|1972l ) prescription to parametrize the ther- 
mohaline mixing as a diffus i on pr ocess. In contrast to 
IDenissenkov fc VandenBer^ (|2003f) . they find that the 
evolution of the surface abundances of carbon, nitrogen, 
and lithium are all affected by overall stellar metallicity, 
while the ^^C/^^C ratio approaches its equilibrium value 
very quickly at all metallicities. 

Although questions of deep mixing rate (e.g.. Smith & 
MarteU 2003) and depth (e.g., Charbonnel et al. 1998) 
have been studied observationally by many authors, the 
results available in the literature can be difficult to syn- 
thesize into a single conclusion. Many studies focus on 
one or two particular clusters (e.g.. Da Costa fc Cot- 
treU 1980, Suntzeff 1981, Trefzger et al. 1983, Lee 
1999), or attempt to correlate deep mixing with other 
clust er properties s uch as horizontal branch morphol- 
ogy (iCavallo fc Naga r 2000), stellar rotat ional velocit y 
([Chaname et al.ll2005[) or cluster ellipticity (|Norrisll987D . 
Individual authors and collaborations develop their own 
analysis tools, and the differences between spectral in- 
dex definitions, model atmospheres, spectral synthesis 
engines, and abundance determination methods produce 
significant systematic differences in different authors' 
abundance scal es, as i s clear from literature-compilation 
studies such as ISmithI ([2002). 

One can construct a phenomenological picture of deep 
mixing from this heterogeneous information, and it goes 
roughly as follows: all stars with mass less than ~ 2M0 
will at some point have their hydrogen-burning shell cross 
the ^-barrier. The /i-barrier is larger than the criti- 
cal /i-gradient for deep mixing, so its destruction per- 
mits deep mixing to begin. Deep mixing occurs con- 
tinuously, and involves all material outside the radius 
where the /i-gradient within the outer H-burning shell 
is critical. The onset of deep mixing happens lower on 
the giant branch for higher-metallicity clusters, because 
their /i-barrier is at smaller radius. However, in higher- 
met allicit;Y_sta£S_th£^ shell is more com- 
pact ([Sweigart fc Mengeilll979| ). so that the radius where 
the /x-gradient is critical is relatively further out in the 
hydrogen-burning shell. This means that the material 
mixed to the surface in higher-metallicity stars is less 
processed than in low-metallicity stars. Various authors 
(e.g., Charbonnel et al. 1998, Cassisi et al. 2002) use 
this relation between metallicity and mixing efficiency to 
study the structure of the hydrogen-burning shell. 

Our goal in this project is to determine the relative 
efficiency of deep mixing across a broad range of metal- 
licity by measuring present-day carbon abundances and 
depletion rates from a homogeneous set of globular clus- 
ter red giants in similar evolutionary ph ases. An earlier 
example of this approach is the study of iBell fc Dickenj 
(1980), who found that [C/Fe] on the upper RGB of M3, 
M13 and NGC 6752 correlated with [Fe/H] metalhcity. 

2. THE DATA SET 

In order to construct a data set that could be used 
to compare CH bandstrengths and [C/Fe] abundances 
across a wide range of metallicity, we chose to obtain 
our data exclusively with the Kast double spectrograph 



3 





. , , , 








..V"^"'"'"'^ NGC 


6712 


LMl 1 : 
: 




'^^^'''"'^ M5 


-39 


: 

: 




M2 


-104 


: 
: 


- / / / 








: / ^/ 
^^^^ ^,-wr*wv--.^X 


NGC 


4147 


IV- 13 : 


, VV / 










M56 


-10 




^ VV / 










WIS 


K77 














Fig. 1. — Selected spectra from our data set. My ~ — 1.5 for all 
stars, [Fe/H] varies from -2.25 (M15 K77) to -1.01 (NGC 6712 
LMll) in roughly equal steps. The vertical bars with widely spaced 
shading mark the bandpasses for the comparison bands of S2{CH), 
and the narrowly spaced shading marks the corresponding science 
band. 

on the Shane 3m telescope at Lick Observatory, devot- 
ing 21 nights between July 2004 and August 2006 to the 
data collection. Using a mirror in place of a dichroic, 
we directed all light to the blue side of the spectro- 
graph, where the 600/4310 (moderate-resolution) grism 
produced a pixel spacing of ~ l.sA/pix and a resolu- 
tion of approximately 5.4A over a wavelength range of 
3400 to 5400A. The detector at the time was a thinned 
1200 X 400-pixel Reticon CCD. Table 1 lists names, po- 
sitions, distance moduli, reddenings, metallicities, and 
number of stars observed for each globular cluster in- 
cluded in the survey. Metallicities, distance moduli and 
reddenings were all taken from the February 2003 revi- 
sion of the online compilation of Harris (1996). Photom- 
etry for the individual globular clusters was taken from 
a combination of original photographic color-magnitude 
diagram work and prope r- motion mem bership s t udies : 
ISandase & Walker (1955) for NGC 4147, Sandage' (l95^ 
and Cudworth (1979b) for^MS^E^ (jl955) and Cudwort: 
(|T979ci) for M 5. lArd (1195511 and : Cudworthl (11979" ^ 
for M13, l Arol (I1955D and 'Harr is eFaL (.1976) for MIO, 
audworth"'|i976bD for M92, Cu dworthI (Il976af) for M15. 



LiUcr (1980) for NGC 6535,|Sandage & Smith' fl966) and 
Cridworth (1988) for N GC 6712,iB arbon (1965) for M56, 
andlHarri^ (|1975D and ICudworth fc Rauscheij (|1987f) for 
M2. 

Targets were observed with a slit width of 1", 1.5", or 
2" depending on the seeing. Ideally, three consecutive 
exposures of each star were taken, with typical exposure 
times of 1800s, to allow for cosmic ray removal. For 17 of 
the 42 stars in our sample, there was only time to obtain 
one or two exposures on a particular night: four of them 
were observed again at a later date, nine have only two 
exposures, and four have only one exposure. In all cases, 
the various exposures were reduced independently, and 
spectra of each individual star were coadded to produce a 
final (unfluxed) 1-d spectrum. The standard stars Feige 
34, Feige 110, BD +28 4211, and BD -^33 2642 were 
observed with a 9" slit aligned with the parallactic angle 
in an effort to capture maximal UV flux for reliable flux 
calibrations. Since Kast sits at Cassegrain focus, spectra 
of the HeHgCd lamp were taken directly following target 



spectra, with the same telescope pointing, slit width and 
dispersive element, to facilitate wavelength calibration 
and account for flexure. 

Data reduction was accomplished with the XIDL 
LowRedux package made available by J.X. Prochaska at 
UC Santa Cruz^. This comprehensive, up dated version 
of the XIDL code used in earlier form in iMartell et al.l 
(|2008bD handles bias subtraction, flat-fielding, cosmic ray 
removal, object identification and extraction, sky sub- 
traction, flexure correction, wavelength calibration, at- 
mospheric correction, coadding and flux calibration. 

Our targets are bright red giants in the range — 1.0 > 
My > —2.0, which generally required 3600-second expo- 
sures with Kast to obtain signal-to-noise ratios per pixel 
just redward of the G band of roughly 150. These stars 
are all significantly brighter than the RGB luminosity 
function bump, meaning that any reasonable deep mix- 
ing rate will have had time to make a measurable decre- 
ment in surface carbon abundance. Figure [1] shows se- 
lected spectra from our sample, with metallicities ranging 
from [Fe/H]= -2.29 (lowest spectrum) to [Fe/H]= -1.01 
(highest spectrum) in roughly equal steps. As [Fe/H] 
rises, there is a clear increase in the strength of the Mg 
b and MgH features near 5 170 A, a slight increase in the 
strength of the broad CN absorption feature at 4215A, 
and a reddening of the overall continuum shape that re- 
duces the apparent depth of the Ca II H& K lines at 3935 
and 3970A and the CN bandhead at 3883A. 

3. ANALYSIS 

We use the inde x 5*2 (CiJ), recently defined in 
IMartell et al.l ()2008a[ ) to be sensitive to carbon abun- 
dance and relatively independent of nitrogen abundance 
over a wide range of metallicity, to quantify the strength 
of the CH G band in all of our combined, flux-calibrated 
spectra. As with most spectroscopic indices, 5*2 (Ci?) is 
measured as the magnitude difference between the inte- 
grated flux in the relevant absorption feature (the "sci- 
ence band" ) and the integrated flux in two nearby rela- 
tively absorptio n-free bands ( t he "co mparison bands"). 
As discussed in IMartell et al.l (|2008a[ ). the science band 
for 82(0 H) runs from 4297A to 4317A, and the compari- 
son bands run from 4212A to 4242 A and 4330 A to 4375 A. 
The bandpasses for S2 (CH) are shown as shaded regions 
in Figure 111 with the more widely spaced lines marking 
the comparison bands, and the more closely spaced shad- 
ing lines marking the science band. 

One-sigma error s on measured values of 5*2 (Cff) were 
determined as in IMartell et al.l ()2008bD : for stars ob- 
served three or more times, as is calculated as the stan- 
dard deviation on the mean of the individual index val- 
ues measured from flux-calibrated, uncombined spectra, 
and for stars observed twice, as — ^ Stars 

observed only once are assumed to have errors in their 
S2{C'H) values equal to the mean value of as, which is 
0.0047. 

To convert bandstrengths to [C/Fe] ab undances, we 
follow a method similar to that described in lMartell et al.l 
()2008bf l. matching measured values of 82(0 H) to values 
derived from synthetic spectra interpolated to match the 
data in My and [Fe/H] (the "model grid"). We select the 

^ available from |http:/ /www.ucolick.org/~xavier/LowRedux/| 



Fig. 2. — Error in calculated [C/Fe] abundance resulting from 
error in the measured index S2{CH), for each star in our sample. 
The abundance error o-jQ/p^j is never larger than 0.12 dex, and is 
almost always smaller than 0.05 dex. 

subset of the model grid with a canonical [N/Fe] value 
of +0.6 and interpolate between model S2{CH) values 
to find a preliminary [C/Fe]. We then interpolate the 
model grid to that preliminary [C/Fe] value, allow [N/Fe] 
to have its full range of possible values, and interpo- 
late betw een model va l ues of the CN bandstrength index 
5'(3839) (jNorris et al.lll98i[ ) to calculate a prehminary 
[N/Fe] abundance. If that preliminary [N/Fe] is more 
than 0.1 dex lower or higher than the canonical value, 
wc interpolate the model grid to match it and repeat 
the S'2(Ci/)-matching process to obtain a final [C/Fe] 
abundance. As a check, we also re-calculate [N/Fe] us- 
ing the final [C/Fe] value, and we find that the differ- 
ence between preliminary and final [C/Fe] and [N/Fe] is 
never large. This is due in large part to the nitrogen- 
insensitivity of S2{CH): for a fixed [C/Fe], [Fe/H] and 
Toff, varying [N/Fe] by large amounts does not change 
82(0 H) significantly. 

The synthetic s pectra employ e d are quite similar 
to those used in iMartell et al.l ()2008bf ): values for 
Teff and logg were t aken from 12-Gyr Yale-Yonsei 
(jPemargue et al.l I2004D isochrones calculated for each 
cluster metallicity, and MARCS model atmospheres 
(|Gustafsson et al.l ll975f ) and the SSG spectral synthe- 
sis program (Bell, Paltoglou, & Trippico 1994 and refer- 
ences therein) were used to generate synthetic spectra. 
For each individual globular cluster metallicity, [C/Fe] 
varies from —1.4 to -1-0.4 in steps of 0.2 dex, and [N/Fe] 
varies from —0.6 to -1-2.0 in steps of 0.2 dex. As in 
IMartell et all (|2008bf) . other variables such as ^^C/^C, 
[0/Fe], and V turb were chosen to be consistent with the 
values used in lBrilev fc Coheiil (|200lD . Synthetic spectra 
were smoothed to a resolution of 5. 4 A and a pixel spacing 
of I.SA to match the data. 

The errors on our [C/Fe] determinations have several 
sources, some of which are readily quantifiable: the er- 
rors in the measured values of ^2 (CH) are generally quite 
small, and propagate into small noise-based [C/Fe] er- 
rors. We determine these noise-based [C/Fe] errors for 
each star in our sample as half the difference between 
[C/Fe] calculated for S2iCH)+as and [C/Fe] calculated 
for 82(0 H)— as- As can be seen in Figurej^l these errors 
are generally smaller than 0.05 dex, and are all smaller 



[Fe/H] 

Fig. 3. — Method-based error A[C/Fe] as a function of [Fe/H], for 
synthetic spectra with -2.31 <[Fc/H]< -0.83, -1.0 <[C/Fe]< 0.0, 
0.0 <[N/Fe]< 1.0, and 2% Poisson noise. Method-based errors for 
individual stars in our sample were assigned as the RMS of all 
points in this figure within ±0.05 o f the [Fe/H] abundanc e of the 
parent cluster. Figure adapted from IMartell et al.l l)2008al ). 

than 0.12 dex. They are also comparable in magnitude to 
the errors introduced by our index-matching method (the 
"method-based error" ) . We calculate a method-based er- 
ror by choosing a random ([Fe/H], [C/Fe], [N/Fe]) point 
within the abundance range spanned by our model grid, 
interpolating to find the 82(0 H) value that corresponds, 
and using that value to re-calculate [C/Fe] assuming a 
canonical [N/Fe] value of -1-0.6. The difference between 
the original input value of [C /Fe] and the calculated value 
is the me thod-based er r or. It is discussed more thor- 
oughly in IMartell et al.l (|2008a[ ). and is a complex func- 
tion of all of the input abundances . Figure [31 which is 
adapted from IMartell et aTl (|2008a| ). shows model-based 
error for a set of synthetic spectra with 2% Poisson- 
distributed noise and abundances in the range typically 
inhabited by globular cluster stars (—2.31 <[Fe/H]< 
-0.83, -1.0 <[C/Fe]< 0.0, 0.0 <[N/Fe]< 1.0), as a 
function of [Fe/H]. The addition of noise to the syn- 
thetic spectra was done to make them resemble the ob- 
served spectra they are compared to, so that the calcu- 
lated method-based error would be an accurate repre- 
sentation of the actual error introduced by our carbon- 
determination method. We assign a method-based error 
for a given star in our survey as the RMS of all points in 
Figure[3]that faU within ±0.05 dex of the [Fe/H] metal- 
licity of its parent cluster. 

We combine the noise-based error and the method- 
based error on calculated [C/Fe] in quadrature, since 
they are independent, and find that the combined er- 
ror ac never exceeds 0.12 dex, and is usually under 0.05 
dex. Table 2 lists cluster name, star name, date observed, 
cluster [Fe/H], My, 82(0 H), as, [C/Fe], and ac for each 
star in our sample. It should be noted that we do not 
include possible errors in our adopted values of [Fe/H] or 
the corresponding T^s and \og(g) values taken from the 
iDemarque et al.l (|2004( l isochrones in ac- 

The matter of systematic offsets to our abundance scale 
resulting from our flux calibration, temperature scale, or 
choice of model atmosphere is not as straightforward to 
assess. We compiled Figure |4| to explore the relationship 
between [C/Fe] values calculated by using a number of 
different G-band indices defined in the literature. Each 



5 



-0.8 -0.6 -0.4 



-0.8 

[C/Fe] froi 



0.6 -0.4 
S,(CH) 



Legend 

)^[C/Fe 
A From 
□ From 
+ From 



] from Set. 

CHfG) 
SfCH) 



-0.8 -0.6 -0.1 



Original data 








Flux — adjusted 








+ 












+ 




















+ 


o 






o ft 

+ + 


+ 


















o - 


ft 




+ 

O 
O 




A 






-6" 




A 


□ 


□ 






□ 












□ 

1 


□ 






Legend 




.1 


. . . . 1 , . 




ft [C/Fe] from s,. 

■' ■ From mcH + Frorr 

n From CH(G) * Frorr 


1 S(CH) 
1 S(4243) 


-0.9 


-O.S 


-0.7 




-0.9 -0.8 


-0.7 



[C/Fe] from S,(CH) 



Fig. 4. — Comparison of [C/Fe] values calculated from various 
CH bandstrengtii indices, for six of tlie eleven clusters in our sam- 
ple. There is generally a vertical offset between abundances derived 
from the various indices, which may be indicative of errors in flux 
calibration. 

vertical column of points represen ts one star, with [C/Fe ] 
calcul ated from the indi ces Sch (iBriley fc Smithlll993D . 
mcH (iSmith et al.lll996n . and CH{G) (|Led[l99^ (which 
were created for studies of red giants in the moderate- 
metallicity globular clusters M13, M3, and M5) plotted 
as stars, triangles, and squares, respectively. Crosses and 
diamonds, respective ly, represent [C/Fe] c alculated from 
the indices S (CH) ("M artell et all l2008br i and S'(4243) 
(jBrilev et al.l fl990). which were designed to be used with 
red giants in the low-metallicity globular clusters M53, 
M55 and NGC 6397. The comparison and science band- 
passes of these indices are given in Table 3. To deter- 
mine these alternate [C/Fe] values, we followed an index- 
matching process identical to the one we used to calculate 
[C/Fe] from 82(0 H). 

Interestingly, in Figure HI there tend to be consistent 
vertical offsets between [C/Fe] values from the various 
indices. Since each index covers slightly different wave- 
length regions, we interpret this to mean that each one 
captures different information from the spectrum. In 
general, indices tuned for low-metallicity stars can have 
wider comparison bands, and do not need to be as care- 
fully defined as indices for high-metallicity clusters, be- 
cause there are absorption features that are important 
to avoid in high-metallicity spectra (for example, the 
CN banhdead at 4215A). This makes the broader lower- 
metallicity indices less affected by noise in the observed 
spectra, but unreliable and often nitrogen-sensitive when 
used on high-metallicity spectra. 

In Figure m the values for [C/Fe] calculated from sch 
are almost always the highest value for any particular 
star. Since sch is defined to have only one comparison 
band, it is the most susceptible of all the indices we mea- 
sure to errors in flux calibration or temperature scale. To 
investigate the possible effects of a mismatch between the 
overall shape of the data and the synthetic spectra, we 
performed a test on our M13 spectra, which have a fairly 
large range in [C/Fe] calculated from the various indices. 
For each of the five M13 stars in our sample, we adjusted 
the flux calibration to match the overall spectral shape of 
an appropriate synthetic spectrum. We fit a parabola to 
fairly absorption-free regions between 4000A and 4500A 



Fig. 5. — Results of the flux-adjustment test: vertical scat- 
ter among [C/Fe] values calculated for M13 stars from various 
CH-sensitive indices is greatly reduced when the overall shape of 
the data continuum is matched to the appropriate synthetic spec- 
trum. The left panel shows [C/Fe] calculated from five G-band 
indices defined in the literature (y-axis) versus [C/Fe] calculated 
from 82(0 H) (x-axis) for the five M13 stars in our data set. The 
right panel shows the results of the same measurements made on 
data adjusted to match the overall shape of abundance-matched 
synthetic spectra. Large circles in the right panel show the [C/Fe] 
values calculated from 82(0 H) and the non-adjusted data spectra, 
and lie very close to the dashed line showing 1:1 correspondence. 

in both the data and interpolated-synthetic spectra, and 
used the ratio of those two curves to adjust the overall 
shape of the data to match the synthetic spectra. The 
left panel of Figure[5]is the same as the JVI13 panel of Fig- 
ure 01 and shows the [C/Fe] abundances calculated from 
the various CH indices. The right panel shows [C/Fe] 
abundances calculated from the same CH indices, mea- 
sured from the fiux-adjusted spectra, and the vertical 
scatter is considerably decreased. The large circles in 
the right panel show the S'2 (Ci? )-derived [C/Fe] values, 
measured from the original (non-fiux-adjusted) spectra, 
and they fall almost exactly on the dashed line showing 
1:1 correspondence. We make two conclusions from this 
exercise: first, that disagreement between [C/Fe] calcu- 
lated from indices that purport to measure the same un- 
derlying abundance can be a sign of flux-calibration or 
temperature-scale errors, and second, that continuum- 
division is likely preferable to flux-calibration, so that 
this potential source of error can be avoided altogether. 

However, the [C/Fe] values calculated from 82(0 H) 
ar e barely chang e d by th e flux adjustment. As discussed 
in iMartell et al.l (|2008aD . 82(0 H) was designed specif- 
ically to be fairly sensitive to carbon abundance and 
insensitive to nitrogen abundance over a wide range in 
metallicity. Therefore we feel that our [C/Fe] values are 
fairly robust, and we choose to use 82ICH) for carbon 
determinations in all of our data, rather than building a 
patchwork of metallicity-tuned indices and correcting for 
estimated systemic offsets. In addition, since our goal is a 
differential measurement of [C /Fe] between stars of vary- 
ing [Fe/H], the absolute zeropoint of our [C/Fe] abun- 
dance scale is not vital to our result. 

In Figure [5] we plot [C/Fe] derived from our 82{CH) 
measurements versus [C/Fe] values taken from the liter- 
ature, for the eight stars in our sample with published 
carbon abundan ces. Carbon abu ndances in JVI3 and M13 
were taken from ISuntzefil ()1981[ ). carbon abundance val- 
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[C/Fe] (this work) 

Fig. 6. — Comparison of our [C/Fe] values calculated from 
82(0 H) with literature ca rbon abundance s. Literature data 
for M3 and M13 are fro m ISuntzefj lUMlh . the MIO data are 
Smitl^_e^^.' (2005"), and the M92 data point is taken 
Bellman et al. (2001). The dashed line shows the line 
which is consisten t wi th the 
BrilevI (|2006l ) and lSmithI 



from 
from 
[C/Fe] 

thiswork 

= [C/Fe]iit + 0.4 

literat ure-compilation studies of ISmith 



(I2OOI) . 



ues for MIQ are taken from [Smith et al.' ('2005'), and the 
single point for JVI92 is from Bellman et al. (2001 ). The 
data compilations of Smith & BrilevI (|200 6) and iSmithI 
([2002) find that Suntzefi (1981) carbon abundances tend 
to fall ~ 0.35 dex below other reported carbon abun- 
dances in 1VI13 and ]VI3. It is therefore encouraging to 
find that the dashed line in Figure [6l which shows the 
relation [C/Fejut =[C/Fe]thiswork — 0.4, traces the data so 
closely. As for the IVIIO and 1VI92 stars, our [C/Fe] values 
are clearly larger than those reported in the literature, 
but with such small crossover samples it is difficult to 
describe this decisively as a systematic offset. 

4. RESULTS 

A functional form for the evolution of carbon abun- 
dance with time may be found by considering the stellar 
envelope as a simple system with a constant mass M and 
mixing rate M . If deep mixing removes carbon from the 
combined convective envelope plus atmosphere at a rate 
—My. Xc,env (where Xc,env is the mass fraction of car- 
bon in the envelope at time t) , and introduces it at a rate 
Al X Xc,hbs (where Xc,hbs is the mass fraction of carbon 
in the hydrogen-burning shell), then the rate of change 
of atmospheric carbon abundance can be written as 



djMXc^env) 

dt 



M{X, 



C,hbs 



X, 



C,env } 



The equilibrium abundance of carbon in the CNO cycle 
is quite small, so this integrates simply to 



X, 



C,env OC CXp 



where M/M can be thought of as the characteristic 
timescale for mixing. Since [C/Fe] is a logarithmic 
measure of carbon abundance, an exponential decline 
in Xc ^nv will result in a line ar decline in [C/Fe] with 
time. ISmith fc IVTartel]! (|2003f ) find a linear decline in 
[C/Fe] with My among bright red giants in the field 
and in the globular clusters 1VI92, NGC 6397, and 1VI3, 
with a slope of 0.22 dex/magnitude. The field-star 



S,(CH) 

Fig. 7. — Derived [C/Fe] abundance versus measured S2{CH) 
bandstrength for all stars in our survey, with total error bars as 
described in the text and symbols denoting three broad [Fe/H] 
bins. Circles represent stars with [Fe/H]< —1.7, triangles are 
for the bin —1.7 <[Fe/H]< —1.4, and squares show stars with 
[Fe/H]> —1.4. Equal [C/Fe] abundances are expressed as different 
CH bandstrengths depending on [Fe/H]. 
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Fig. 8. — Derived [C/Fe] abundance with error bars, as in 
Figure [7] versus [Fe/H]. Here symbols denote three broad My 
bins. Circles show stars with My > —1.3, triangles are for the bin 
— 1.3 > My > —1.6, and squares are for stars with My < —1.6. 
Although the faintest and brightest bins separate fairly well into 
high and low carbon abundances, respectively, it is unclear from 
this plot whether low-metallicity stars experience more rapid car- 
bon depletion than high-metallicity stars. 



study of iGratton et al.l ()2000l ) finds a linear trend in 
[C/Fe ] with log(L/L0). According to iDemargue et all 
(|2004f ) evolutionary tracks, the \og{L / LQ)-t relationship 
is not exactly linear for post-bump RGB stars (1.8 < 
log(i/LQ) < 2.6), but it is sufhciently close that a linear 
[C/Fe] - My or [C/Fe] - log(L/L0) relation suggests a 
fairly linear [C/Fe] - t relation as well. 

Figure [7] shows calculated [C/Fe] versus measured 
82(0 H) for our sample, with symbols indicating three 
broad metallicity bins. Circles are stars from clus- 
ters with [Fe/H]< —1.7, triangles represent stars with 
— 1.7 <[Fe/H]< —1.4, and squares are for stars with 
[Fe/H]> —1.4. It can be clearly seen that equal carbon 
abundances are expressed as lower CH bandstrengths at 
lower metallicity. The error bars shown are equal to ±ac- 
Figure [H] shows calculated [C /Fe] with error bars as a 
function of [Fe/H], and here the symbols represent three 
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Fig. 9. — The relation between the absolute My magnitude 
of the "bump" in the red giant branch luminosity function and 
metallicit y. Large circles are data from 11 globular clusters ob- 
served in lFusi Pecci et al.l ||T990), and small circles are the values 
for the globular clusters in our sample, interpolated linearly from 
the IFusi Pecci et al.l | |T990|) data. 

broad luminosity bins. These bins have the same sense 
as in Figure [71 with the faintest stars represented by cir- 
cles and the brightest stars plotted as squares. Although 
the faintest and brightest stars occupy mostly separate 
regions in this figure, there is not a clear trend between 
[C/Fe] and [Fe/H]. However, the lack of a clear correla- 
tion in Figure [5] does not necessarily mean that there is 
no dependence of deep mixing rate on stellar metallicity. 

The question of how much time each observed star 
has spent in the deep mixing phase is the key to 
interpreting our data: we selected targets with sim- 
ilar absolute magnitudes, but it is well-known ob- 
servationally jF usi Pec ci et al.l |199C1[ ) and theoretically 
(jCassisi fc Sala ris 1997) that the RGB bump occurs at 
lower luminosity in more metal-rich globular clusters. In 
addition, low-metallicity red giants evolve more quickly 
than high-metallicity red giants of the same mass (e.g., 
Riello et al. 2003). As a result, the My = -1.5 red gi- 
ants in the metal-rich end of our sample have been mixing 
for far longer than their metal-poor counterparts. In or- 
der to study deep mixing efficiency with our data set, we 
must convert our [C/Fe] values, which are a function of 
both [Fe/H] and time, to some time- invariant quantity. 
We choose to do this by converting present-day [C/Fe] 
abundances to a carbon depletion rate per Gyr. This is 
a two-step process: we must first calculate the absolute 
magnitude of the RGB bump, M^™'^ , for each individ- 
ual cluster, and then for each star we must convert the 
V magnitude height above the bump, AM^"™^, into a 
time since the onset of deep mixing. 

We use the observational study of iFusi Pecci et al.l 
(|1990f ) to calculate M^"™^ for each cluster in our sur- 
vey. As can be seen in Figure [9l there is a nearly- 
linear relationship between observed My"™^ (plotted as 

large open circles) and [Fe/H]. We calculate My"™^ for 
the eleven clusters i n our sample by a simple linear in- 
terpolation from the iFusi Pecci et al.l (fT99Q) data. The 
conversion from AM,''""^ to At''""^P is based on Yale- 
Yonsei (jPemarque et al.ll2004[ ) isochrones and evolution- 
ary tracks. For each cluster metallicity, we created an 
evolutionary track for a star with appropriate metal- 
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Fig. 10.— Rate of evolution AMy/At versus [Fc/H] for the 
eleven globular clusters in our sample, calculated using Yale-Yonsei 
(Dcmarquc ct al. 2004) isochrones and evolutionary tracks in the 

range M^"^^ > My > —1.5. The solid curve is the best-fit poly- 
nomial, and the dashed curve is the best-fit exponential. 

licity and a mass taken from the My — —1.5 point 
in a metallicity-matched 12 Gyr YY isochrone. We 
then chose points in the evolutionary track near My"™^ 
and My = —1.5, and we calculate AMy/At based on 
those points. Figure [TOl shows our calculated values for 
AAfy/Ai for all clusters in our sample, along with a 
best-fit polynomial (solid curve) and a best-fit exponen- 
tial (dashed curve), which are nearly identical. 

Combining these two steps, we calculate the carbon de- 
pletion rate A[C/Fe]/At as [C /Fe] / AMy"'^ x AMy/At, 
assuming that the initial [C/Fe] abundance is solar, and 
that deep mixing begins at the RGB bump for all stars. 
Figure [Til shows A [C/Fe] /At versus [Fe/H] for all stars 
in our sample, and there is a clear downward trend in 
carbon depletion rate as metallicity increases. We see no 
break or corner in the relation between A[C/Fe]/Ai and 
[Fe/H], indicating that mixing is never completely pro- 
hibited in our sample. This may happen at metallicities 
above [Fe/H]= —1.0, if the hydrogen-burning shell is so 
compressed that the critical /^-gradient is reached out- 
side it. Our r esult is in agreement with the theoretical 
predictions of [Charbonnel fc ZahnI (|2007| ). and confirms 
that the process of deep mixing is less effective in rela- 
tively high-metallicity red giants. 

5. CONCLUSIONS AND FUTURE WORK 

In summary, the results of this work show that within 
the [Fe/H] range —2.2 to —1.0 dex the rate of deep mix- 
ing varies with metallicity. We measure CH bandstrength 
using the index 82(0 H), which w as designed to be vali d 
across a broad range in [Fe/H] ([Martell et al.l[2008ah . 
Carbon abundances are determined by matching CH 
bandstrengths measured from the data to bandstrengths 
measured from specifically-designed grids of SSG syn- 
theti c spectra. Under the assumption ([Charbonnel et all 
[1998( 1 that deep mixing begins at the RGB luminosity 
function bump, we establish the carbon depletion rate 
for a given star as the change in its [C /Fe] from an initial 
solar value divided by the time elapsed since the onset 
of mixing. Since the RGB luminosity function bump is 
fainter in high-metallicity globular clusters than in lower- 
metallicity globular clusters, and higher- metallicity red 
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models of deep mixing. Nevertheless, our fundamental 
result, that carbon depletion proceeds more quickly in 
low-metallicity globular cluster red giants than in their 
high-metallicity counterparts, is robust, and provides a 
clear affirmation of present theoretical models. 
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Fig. 11.— Carbon depletion rate A[C/Fe]/At versus [Fe/H] for 
all of the stars in our sample. There is a clear downward trend with 
risin g metallicity, a result predicted by the deep mixing models of, 
e.g.. lCharbonnel fc Zahiil ((2007h . 

giants evolve more slowly than lower metallicity red gi- 
ants, highcr-metallicity stars at a given absolute magni- 
tude qualitatively ought to have spent more time mix- 
ing than their lower-metallicity counterparts. In or- 
der to quantify the time spent mixing as a function of 
metallicity and absolute V magnitude, we interpolate 
the lFusi Pecci et"al] (|1990f ) observational data on M^"™^ 
as a function of metallicity to the mctallicities of the 
clusters in our study, then use m etallicity-appropriate 
Yale-Yonsei ()Demaraue et al.l2004f) isochrones to convert 
^j^bump g^^j^ individual star into At''"™P. As can 
be seen in Figure (TOl the lower-metallicity red giants do 
evolve more quickly, and therefore spend less time mixing 
than the higher-metallicity red giants in this study. How- 
ever, dividing [C/Fe] by Ai''"™^ for each individual star 
produces the result that the deep mixing rate is roughly 
twice as large at [Fe/H]= —2.0 than at [Fe/H]= —1.0, as 
can be seen in Figure [TTl 

This present analysis includes some assumptions worth 
noting, since they may need to be more carefully ex- 
amined if this result is to provide constraints for theo- 
retical models of deep mixing. To express Xc,env as a 
straightforward exponential in time, we hold the mass of 
the stellar envelope constant, though it shrinks contin- 
uously as the hydrogen-burning sh ell proceeds outward 
(|Denissenkov fc VandenBerd 120031 ) . We also hold the 
mixing rate M constant, though it is controlled by the 
structure of the hydrogen-burning shell, and may well 
evolve. In addition, we implicitly assume in construct- 
ing Figure [TT] that all stars have equal (and solar) ini- 
tial abundances of carbon. A primordial depletion of 
~ 0.3 dex in a subset of our sample would thus be mis- 
interpreted as an artificially high mixing rate in those 
stars. A study of carbon abundances in fainter post- 
bump red giants, or pre-bump giants, in the globular 
clusters included in this study, would allow for direct cal- 
culation of A [C/Fe] /At without assumptions about the 
initial carbon abundance. In addition, a study of car- 
bon depletion rates in high-metallicity red giants, either 
in old open clusters or in high-metallicity globular clus- 
ters, would allow an estimation of the maximum metallic- 
ity at which deep mixing still operates. That maximum 
metallicity would be helpful for constraining the various 



TABLE 1 



Cluster ID RA (J2000) S (J2000) (m - M)v E{B - V) [Fe/H] N^bs 



NGC 


4147 




12 


10 


06.2 


+18 


32 


31 


16.48 


0.02 


-1.83 


4 


NGC 


5727 


(M3) 


13 


42 


11.2 


-28 


22 


32 


15.08 


0.01 


-1.39 


3 


NGC 


5904 


(M5) 


15 


18 


33.8 


+02 


04 


58 


14.46 


0.03 


-1.29 


3 


NGC 


6205 


(M13) 


16 


41 


41.5 


+36 


27 


37 


14.45 


0.02 


-1.54 


5 


NGC 


6254 


(MIO) 


16 


57 


08.9 


-04 


05 


58 


14.08 


0.28 


-1.52 


8 


NGC 


6341 


(M92) 


17 


17 


07.3 


+43 


08 


11 


14.64 


0.02 


-2.29 


2 


NGC 


6535 




18 


03 


50.7 


-00 


17 


49 


15.22 


0.34 


-1.80 


2 


NGC 


6712 




18 


53 


04.3 


-08 


42 


22 


15.60 


0.45 


-1.01 


3 


NGC 


6779 


{M56) 


19 


16 


35.5 


+.30 


11 


05 


15.65 


0.20 


-1.94 


5 


NGC 


7078 


{M15) 


21 


29 


58.3 


+ 12 


10 


01 


15.23 


0.10 


-2.25 


1 


NGC 


7089 


(M2) 


21 


33 


29.3 


-00 


49 


23 


15.49 


0.06 


-1.62 
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TABLE 2 



— 

Cluster ID 


— 

Star ID 


Date obs. 


My 


[Fe/H] 


S2{CH) 






[C/Fc] 






MIO 


11-105 


2006- 


■06-03 


-1.23 


-1, 


.52 


1. 


.753 


0, 


,0114 


-0, 


,625 





,0549 


MIO 


III- 73 


2006- 


■06-03 


-1.22 


-1 


.52 


1. 


.784 


0, 


,0011 


-0, 


,507 





,0327 


MIO 


III-85 


2006- 


■06-03 


-1.38 


-1 


.52 


1. 


.784 


0, 


,5153 


-0, 


,544 





,0325 


MIO 


IV-30 


2006- 


■06-03 


-1.28 


-1 


.52 


1 


.803 


0, 


,0052 


-0, 


,456 





,0381 


MIO 


1-15 


2004- 


■07-12 


-1.22 


-1, 


.52 


1, 


.799 


0, 


,0047 


-0, 


,461 





,0365 


MIO 


III-93 


2004- 


■07-12 


-1.18 


-1, 


.52 


1, 


.802 


0, 


,0047 


-0, 


,448 





,0370 


MIO 


III-97 


2004- 


■07-12 


-1.41 


-1, 


.52 


1, 


.724 


0, 


,0047 


-0, 


,764 


0, 


,0375 


MIO 


1-12 


2005- 


■07-12 


-1.14 


-1, 


.52 


1, 


.819 


0, 


,0102 


0, 


,375 


0, 


,0513 


M13 


IV-53 


2005- 


■04-17 


-1.77 


-1 


.54 


1 


.747 


0, 


,0019 


-0, 


,729 


0, 


,0336 


M13 


11-33 


2006- 


■06-01 


-1.78 


-1, 


.54 


1, 


.764 


0, 


,0021 


-0, 


,647 


0, 


,0337 


M13 


III-52 


2006- 


■06-01 


-1.78 


-1, 


.54 


1, 


.735 


0, 


,2635 


-0, 


,763 


0, 


.0325 


M13 


11-57 


2006- 


■06-02 


-1.75 


-1, 


.54 


1, 


.703 


0, 


,0023 


-0, 


,910 


0, 


.0342 


M13 


III-18 


2006- 


■06-02 


-1.68 


-1, 


.54 


1, 


.698 


0, 


,0020 


-0, 


,914 


0, 


.0338 


M15 


K77 


2005- 


■09-06 


-1.53 


-2, 


.25 


1, 


.584 


0, 


,0016 


-0, 


,718 


0, 


.0386 


M2 


1-103 


2005- 


■07-13 


-1.92 


-1, 


.62 


1, 


.775 


0, 


,0071 


-0, 


,582 


0, 


.0447 


M2 


1-104 


2005- 


■07-13 


-1.52 


-1, 


.62 


1, 


.792 


0, 


,0042 


-0, 


,486 


0, 


.0359 


M2 


1-298 


2005- 


■09-07 


-1.46 


-1, 


.62 


1, 


.755 


0, 


,0037 


-0, 


,601 


0, 


,0351 


M2 


1-190 


2005- 


■09-08 


-1.75 


-1, 


.62 


1, 


.698 


0, 


,0024 


-0, 


,888 


0, 


,0340 


M2 


11-60 


2005- 


■09-09 


-1.69 


-1, 


.62 


1, 


.772 


0, 


,0034 


-0, 


,570 


0, 


,0348 


M2 


11-71 


2005- 


■09-09 


-1.64 


-1, 


.62 


1, 


.771 


0, 


,0042 


-0, 


,570 


0, 


,0361 


M3 


BC 


2006- 


■05-31 


-1.24 


-1, 


.39 


1, 


.825 


0, 


,0010 


-0, 


,423 


0, 


,0333 


M3 


1-46 


2006- 


■05-31 


-1.26 


-1, 


.39 


1, 


.809 


0, 


,0050 


-0, 


,488 


0, 


,0389 


M3 


V-80 


2006- 


■06-02 


-1.61 


-1, 


.39 


1, 


.763 


0, 


,0031 


-0, 


,702 


0, 


,0361 


M5 


L39 


2006- 


■06-01 


-1.39 


-1, 


.29 


1, 


.838 


0, 


,0041 


-0, 


,342 


0, 


,0410 


M5 


IV-34 


2006- 


■06-02 


-1.41 


-1, 


.29 


1, 


.791 


0, 


,0045 


-0, 


,584 


0, 


,0418 


M5 


IV-49 


2006- 


■06-03 


-1.32 


-1, 


.29 


1, 


.713 


0, 


,0031 


-0, 


,917 


0, 


,0384 


M56 


I-IO 


2006- 


■08-30 


-1.53 


-1, 


.94 


1, 


.682 


0, 


,0022 


-0, 


,610 


0, 


,0304 


M56 


E-48 


2006- 


■08-31 


-1.95 


-1, 


.94 


1, 


.694 


0, 


,0031 


-0, 


,678 


0, 


.0315 


M56 


1-141 


2005- 


■09-06 


-1.22 


-1, 


.94 


1, 


.733 


0, 


,0025 


-0, 


,343 


0, 


.0301 


M56 


E-22 


2005- 


■09-07 


-1.78 


-1, 


.94 


1, 


.701 


0. 


,0035 


-0, 


,597 


0, 


.0318 


M56 


1-66 


2005- 


■09-07 


-1.78 


-1, 


.94 


1, 


.653 


0, 


,0015 


-0, 


,798 


0, 


.0297 


M92 


11-70 


2006- 


■05-31 


-1.40 


-2, 


.29 


1, 


.621 


0, 


,0036 


-0, 


,398 


0, 


.0420 


M92 


IV-94 


2006- 


■05-31 


-1.44 


-2, 


.29 


1, 


.549 


0, 


,0047 


-0, 


,823 


0, 


,0563 


NGC 4147 


11-14 


2005- 


■02-01 


-1.150 


-1, 


.83 


1, 


.697 


0, 


,0123 


-0, 


,580 





.0539 


NGC 4147 


11-30 


2005- 


■02-01 


-1.760 


-1, 


.83 


1, 


.754 


0, 


,0027 


-0, 


,519 





.0313 


NGC 4147 


11-45 


2005- 


■02-02 


-1.960 


-1, 


.83 


1, 


.668 


0, 


,0093 


-0, 


,913 





.0516 


NGC 4147 


IV-13 


2005- 


■02-02 


-1.390 


-1, 


.83 


1, 


.713 


0, 


,0066 


-0, 


,585 





.0381 


NGC 6535 


13 


2004- 


■07-13 


-1.71 


-1, 


.80 


1, 


.738 


0, 


,2743 


-0, 


,583 


0, 


,0308 


NGC 6535 


19 


2004- 


■07-13 


-0.990 


-1, 


.80 


1, 


.781 


0, 


,0068 


-0, 


,288 


0, 


,0379 


NGC 6712 


KC564 


2005- 


■07-13 


-1.430 


-1, 


.01 


1, 


.850 


0, 


,0147 


-0, 


,335 





.1187 


NGC 6712 


LMll 


2005- 


■07-13 


-1.200 


-1, 


.01 


1, 


.838 


0, 


,0139 


-0, 


,397 


0, 


,0974 


NGC 6712 


B66 


2006- 


■08-31 


-1.57 


-1, 


.01 


1, 


.808 


0, 


,0110 


-0, 


,676 


0, 


,1069 
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TABLE 3 



Index 


Blue comparison band (A) 


Science band (A) 


Red comparison band (A) 


Reference 


SCH 

mcH 

CH(G) 

S{CH) 

5(4243) 

S2{CH) 


4220-4280 
4080-4130 
4230-4260 
4050-4100 

4212-4242 


4280-4320 
4270-4320 
4270-4320 
4280-4320 
4290-4318 
4297-4317 


4420-4470 
4390-4420 
4330-4350 
4314-4322 
4330-4375 


Brilcv & Smith ( 1993) 
Smith et ah (1996) 
Lee (1999) 

MartcU et al (2008b') 
Brilev et aTTiggO) 
Martell et al. (S'OOSa) 
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